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ABSTRACT 

We develop a method for estimating the shear power spectra from weak lensing observations and 
test it on simulated data. Our method describes the shear field in terms of angular power spectra and 
cross correlation of the two shear modes which differ under parity transformations. Two of the three 
power spectra can be used to monitor unknown sources of noise in the data. The power spectra are 
decomposed in a model independent manner in terms of "band-powers" which are then extracted from 
the data using a quadratic estimator to find the maximum of the likelihood and its local curvature (for 
error estimates) . We test the method against simulated data from Gaussian realizations and cosmological 
A-body simulations. In the Gaussian case, the mean bandpowers and their covariance are well recovered 
even for irregular (or sparsely) sampled fields. The mild non-Gaussianity of the A-body realizations 
causes a slight underestimation of the errors that becomes negligible for scales much larger than several 
arcminutes and does not bias the recovered band powers. 

Subject headings: cosmology: theory - gravitational lensing - large-scale structure of universe 



> 

in 
m 
O 

o : 
o ■ 

6 ■ 
a: 

«a . 



1. INTRODUCTION 



Weak lensing of background galaxies by foreground 
l arge-scale structure has now been convincingly detected 
( [Bacon. Rcfrcgicr fc Ellis 2000; Kaiser, Wilson fc Luppinc) 



2000T 



Waerbeke et al 



IVit.tman et al.| 2000) 



206TJ^ 

and is rapidly becoming a valuable tool for studying the 
distribution and clustering evolution of dark matter in the 
universe. Lensing produces a correlated distortion of the 
ellipticities of background galaxies, at the percent level, 
which can be used to measure a two-dimensional projec- 
tion of the mass distribution. 

While there are numerous observables that can be de- 
fined from the shear maps produced by such surveys, one 
of the most important is the angular power spectrum. 
The angular power spectrum contains valuable informa- 
tion on cosmological paramet ers that compl ement s other 



astro physical measure ments ( [Iain fc Seljak 1997; Kaiser 



1998; Hu & Tegmark 1999). In this respect weak lensing 



is very similar to anisotropies in the CMB and in fact much 
of the theory and data analysis is very analogous as well. 

In this paper, we suggest a technique for the presen- 
tation and interpretation of weak lensing data which is 
commonly employed in CMB analysis: the use of "band- 
powers" extracted from the data by an iterated quadratic 
estimator of the maximum likelihood solution. This tech- 
nique has the advantage of automatically taking into ac- 
count irregular survey geometries and varying sampling 
densities. It provides an optimal estimate of the power 
spectrum in the Gaussian regime and makes efficient use 
of all of the data on the relevant angular scales. Error esti- 
mates that include the sampling and noise variance of the 
survey are also automatically produced by this method. 
Similar claims cannot be made for methods based on corre- 
lation functions or simple Fourier Transforms of the data. 

Throughout our focus will be on weak lensing by large- 
scale structure, specifically on angular scales larger than 
several arcminutes (or wavenumbers I <^ 10 3 ). These are 
the scales which will be probed by ongoing wide field sur- 



veys, e.g. the Deep Lens Survey 1 , surveys with the VLT 
and the Hawaii/If A lens survey. On these angular scales 
the signal reduces simply to a p rojection of the density 



contrast along t he lin e of si ght ( Blandford et al] 1991; 
Miralda-Escude 1991; Kaiser 1992), which is well approx- 
imated by a Gaussian probability distribution. In this 
limit the angular power spectrum encodes all of the rel- 
evant information about the field, and in particular can 
be predicted from the 3D power spectrum of the density 
field using Limber's equation. On subarcminute scales, the 
field becomes substantially non-Gaussian, lens-lens cou- 
pling and perturbations to the photon trajectory become 
increasingly important and the weak lensing approxima- 
tions break down. On these angular scales the signal is 
dominated by individual objects along the line of sight 
(e.g. clusters of galaxies) and a correlation function or 
power spectrum analysis becomes less useful. 

The outline of this paper is as follows. In §0 we describe 
the shear signal covariance matrix and its relationship to 
the three fundamental shear power spectra. We then de- 
velop (^) and test (§[!]) the iterated likelihood method for 
band power extraction using simulated data. We conclude 
in §| 

2. WEAK LENSING 

In this section we briefly review the theory of weak grav- 
itational lensing to establish our notation and conventions. 
Throughout we shall focus on the 2-point function of gravi- 
tational shear, since this shall determine the angular power 
spectrum in §||. 

The gravitational deflection of light induces a mapping 
between the 2D source plane (S) and the image plane (I). 
The deformation so induced can be written 



Sxf = A 



■At 1 



(1) 



where 5x is the separation vector between points on the 
respective planes. In the we ak lensin g limi t, the deforma- 
tion can be decomposed as (Mellier 1999; Bartelmann & 

*|http :/ /dls. bell-labs. com 
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Fig. 1. — The fundamental shear modes and their cross correlation, (a) A pure e-field obtained from a convergence map from White & Hu 



(2000). (b) A pure /3-field obtained by a rotation of the shears by 7r/4. (c) A correlated mixture of e and with C" = Cj 
by a rotation of the shears by tt/8. 
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Schneider 2000) 

Aij = (1 - K)5ij - 71CT3 



72CT1 



(2) 



where the <Ji are the 2x2 Pauli matrices, n <C 1 is the 
convergence and 7 Q ^ 1 is the shear. If a galaxy has 
(weighted) second moments M s then the image will have 



M 1 = A- 1 ■ M s ■ A- 1 



(3) 



The ellipticities are usually defined in terms of the sec- 
ond moments of the light distribution, corrected for instru- 
mental and observational effects, and in the weak lensing 
regime, Eq. (0) simplifies dramatically such that the ob- 
served ellipticity of a galaxy is linearly related to the shear. 
The proportionality constant depends on the definition of 
the ellipticity; we take 



( e ) = 7 , 



(4) 



but n ote that 27 i s often found in the literature ( Bartel-| 
mann fc Schneider 2000). The result is that e defines a 



(noisy) estimate of the local shear field at A. 

Now consider an observation of a given area of the sky. 
The observed field yields an estimate of the ellipticities 
ej and positions A^ of a set of galaxies binned into pixels 
i = 1, • • • -/Vpi x . In a Cartesian coordinate system on the 
sky the two components of the shear field, 71(A) and 72(A), 
transform as a spin-2 field. The Fourier decomposition is 



71(A) ± 172(A) 



W(l)[e(l)±i0(l)}e 



±2itpi e iln 



(5) 

where <pi is the angle between 1 and the x-axis and W(l) 
is the Fourier transform of the pixel window function. For 
square pixels of side a in radians 



w ( l ) =3o ( TrCOSi/?; I j 



la 



■ sin ipi 



(6) 



where jo(x) = sin(cc)/a; is the 0th order spherical Bessel 
function. Note that for long wavelengths the pixelization 
is irrelevant and the window goes to unity. 

We shall be interested in the power spectrum or correla- 
tion function of the shear field. The two point correlations 



in the shear are determined by the three shear power spec- 
tra 



(e(l)e(l')) 


= (2tt 


) 2 <5(1 


- i')Q ee 


</W(0> 


= (2tt 


) 2 <5(1 


- l')Cf 


«1)/3(1')> 


= (2tt 


) 2 <5(1 


- \')Cf 



(7) 

For the shear generated by weak lensing Cf e = C, reK , Cf 3 = 



and Ci 



e/3 



0. For shot noise Cf e = C, 



00 



Systematic 



errors can in principle generate any of the power spectra. 

Since a 45° degree rotation of the shears takes e — > /3, 
it converts the lensing signal to a C 00 = Cf K spectrum 
with C" — Cf = 0. A more general rotation leaves 
a signal in both C" and C 00 but also correlates them 
as (Cf* 2 



-,00 



For the shot noise, the relation 



a 



I I - °/ U l 

1 = C 00 is invariant under rotations. These rotations 
also allow one to visualize the pattern implied by each 
spectra (see Fig. |l|). In particular, the (3 component pos- 
sesses a "handedness" ; formally the two are distinguished 
by their transformation under parity. 
By direct substitution 



(7i(Ai)7i(Aj)) 



f—( 

J (2tt) 2 V 



Cf cos 2 2(pi + Cf sin 2 2ip, 



d 2 l 
(2 

Cf sin4 w )w 2 (l)e iI - ( "*- ii ^ 



(72(Ai)72(Aj)) = 



(7i(Ai)72(Aj)) 



d 2 l 
(2^ 



Cf sin 2 2ipi 



Cf cos 2 2(pi 



+Cf sin4^)^ 2 (l)e ll (li ' 



d 2 l rl 



(2nf 



C 00 ) sin4 W 



+C ; e/3 cos4 w j W a (l)e fl -(**-*^ . (8) 

For a coordinate system that is oriented so that n.; — A j \\ 
x and pixel separations that are small compared with the 
coherence scale of the field, the cosmological ee signal gen- 
erates (7171) > 0. For shot noise (7272) ~ (71 7i) and 
for either (7172) ~ 0. These are the tests suggested by 
Miralda-Escudel (1991). 



3 



o.i - 



10-3 = 




Fig. 2. — Low order moments of the pixel window function 
squared. The first two moments sufficiently approximate the pixel 
window function to I ~ 2£ p j x = An/cr where a is the width of the 
pixel in radians. 



These relations (||) allow us to define the lensing signal 
correlation matrix 

Cg )(o6 ) = (7«(fi*)7b(%)> ■ (9) 
The correlation matrix may be simplified by recalling that 

oo 

exp (il • 6) = J Q (W) + 2 i m Jm{lO) cos(m(<^ - (j>)) , 

( 10 ) 

where (0, 0) define the magnitude and orientation of the 
separation vector — iij. Furthermore, the window func- 
tion can be similarly decomposed as 



W 2 (l) = J2w n (l)cos(nip l ). 



(11) 



For square pixels, the n = 1,2,3 moments vanish due to 
symmetry and it is sufficient to retain the isotropic and 
n = 4 quadrupole contributions (see Fig. ||). The integral 
over ifi can now be performed analytically leaving 



sig _ f ldl^-y x 
J x 



1 



(12) 



where X takes on the values ee, (3(3 and e(3 

r e = 



Jo "I" C4J4 S4J4 



J/3/3 = 



Jo — C4J4 — S4J4 
— 54 J4 Jo H - C4 J4 

— 2S4J4 2C4J4 

2c4</4 254^/4 



(13) 



and 



1/3/3 _ 



Q 



e/3 



J + 2c 4 J 4 + C 8 J 8 S 8 J 8 

S 8 J 8 — Jo + 2C4 J4 — C 8 J 8 

—Jo + 2C4 J4 — C 8 Js —sgJg 

—sgjg, Jo + 2c4 J4 — eg Js 

-2s 8 J 8 2J + 2c 8 J 8 
2J + 2c 8 J 8 2s 8 Js 



(14) 



Here we have used the shorthand notation c n = cos(n0) 
and s n = sin(n</>) and the (suppressed) argument of the 
Bessel function in each case is 16. Thus for a fiat band- 
power we need only evaluate 



C {n) {0) 



leB 



dl 
21 



J n (W). 



(15) 



for n = 0, 4, 8. 



3. SHEAR LIKELIHOOD 



We wish to estimate the (angular) power spectra of the 
shear field from the observed image ellipticities by means 
of a maximum likelihood technique. This ensures that, 
under the stated assumptions, we make optimal use of the 
data and correctly handle any irregular survey geometry 
which may affect the correlations on large angular scales. 
Non-uniform or correlated noise can also be efficiently han- 
dled by this formalism. 

First we decide to parameterize the underlying power 
spectra with a set of parameters p a where a = 1, . . . , N p . 
These parameters could be cosmological, or describe a par- 
ticular model of noise or systematic errors. In this paper 
we shall mainly be interested in the case where the p a 
are "bandpowers" , i.e. we shall approximate the angular 
power spectra as piecewise constant with p a the value of 
1(1 + \)Ci/2tt in band B a . So long as the constant sec- 
tions are narrower than any feature in the power spectrum 
which we wish to reproduce, the sharp steps in power will 
not produce any ill effects, and at the same time the pa- 
rameterization is model independent. 

Such a specification is not sufficient to describe the most 
general likelihood function of the observations given the 
p a , however if the shear field is Gaussian it is. We shall as- 
sume that on sufficiently large angular scales, which are of 
interest here, the field is sufficiently Gaussian that the es- 
timate of the 2-point function so derived is not seriously in 
error. Once the field becomes significantly non-Gaussian 
the utility of a power spectrum estimate become suspect. 
We shall return to this point in detail in the next section. 

Consider the data as a 2N component vector 



d = {7i(ni),7 2 (ni); . . . ; 7i(njv), j2(n N )} , 
then the likelihood is simply 



(16) 



£( ^ )= (2^1^)1^ ^- 1 2 dT °- 1 ^- ^ 

Here the correlation matrix C is a sum of two terms. The 
first is the cosmological signal, Eq. (^|), and the second 
is the noise. Each galaxy has an rms intrinsic ellipticity 
per component 7 int which we assume is uncorrelated with 
the underlying shear field taken to be constant across the 
galaxy. The ellipticity is thus a noisy estimator of the 
shear field at its position, and the noise matrix is 



y^noisc 
°(»i)(ab) 



Ni 



-SijS a b 



(18) 



where JVj is the number of galaxies in pixel i. In more gen- 
erality, the noise matrix can include observational errors 
on the ellipticities which could be correlated from galaxy 



4 



to galaxy. The likelihood formalism can efficiently handle 
a general noise matrix. 

One then maximizes the likelihood as a function of the 
model parameters p a . We accomplish this maximization 
iteratively by using the Newton-Raphson me thod to find 



the root of dC/dp a = (Bond, Jaffc & Knox 1998; Seljak 
1998). Specifically, from an initial set p a one makes an 
improved estimate p a = p a + 5p a where 



6 Pa = X^liF- 1 )^ (dd T - C)(C- 1 C^C- 1 ) 



and where the Fisher matrix is 



Fa/3 — 2 tr (^ 



(19) 
(20) 



and A < 1 is a parameter that is set to reduce the step size 
when the full value causes large jumps in parameter space. 
Note that in maximizing this function for the bandpowers 
several matrices, and their derivatives, must be computed 
and inverted. Using Eqs. (12-15|), we first compute the 
derivative matrices C a for p a = 1. The full correlation 

function is then the sum ^ Q p Q C !Q plus the noise term. 

With an appropriate choice of p a and A this method 
typically converges rapidly to the desired solution. Gen- 
eral rules of thumb for these quantities include that the 
bands should be at least twice as wide as /fi c id = 27r/0 max 
where 6* max is the largest pixel separation. This reduces 
the covariance between the bands that is due to the sur- 
vey geometry. The lowest band should include ee modes 
below ifi e id to absorb any d.c. offsets in the data. The 
highest band should include modes above = 27r/# p i X 
since the window function of square pixels has a long tail 
to high multipoles. 

The choice of A to achieve rapid convergence is more 
of an art than a science. We have found empirically that 
it is sufficient to choose A so that at each iteration no 
ee parameter in the range Zg e id < I < ^pix changes up or 
down by a factor of more than 3. However, for bands that 
are noise dominated there are cases where the maximum 
likelihood solution desires negative signal power. For this 
reason, we also impose a minimum A = 0.01 so that the 
power in a band can be reduced below zero; we do require 
that the total power in the signal and noise be positive 
by resetting negative values to a small positive number at 
the start of each iteration. On top of these criteria, we 
halve A whenever 5p's from consecutive iterations cancel 
to » 20%. Since there is in principle no or little signal 
in the (3(3, e/3 and end ee bands we stop iterating when 
we reach better than 5% convergence in all the remaining 
bands. The end bands should then be dropped when using 
the bandpowers for cosmological constraints. 

Finally it remains to estimate the errors on our maxi- 
mum likelihood power spectrum. Ideally one would esti- 
mate the 2-point function of the parameters p a by Monte- 
Carlo integration 



C0V Q/3 = (Sp a Spp) = J dp Sp a Spf3 C(p a ) (21) 

over the likelihood function. However the likelihood eval- 
uation is sufficiently slow that this method impractical. 



Instead we use the Fisher matrix. If the likelihood is suf- 
ficiently Gaussian (in the parameters) around the max- 
imum, one can estimate the covariance matrix from the 
curvature matrix, 



F a = -[ln£(p a )],a/3 

= tr[(dd T - QiC^C^C^C^C- 1 



(22) 



-\c- l C tCt pC- 1 )] + itr(C , - 1 C, a C- 1 C, /9 ) 

or, assuming that that the maximum likelihood model is 
correct, from its expectation value, the Fisher matrix, 



Cov Q/3 w (F r ) a p . 



(23) 



Since the Fisher matrix is automatically calculated in the 
iteration process, error estimates come at no additional 
computational expense. 

4. TESTS 

There are several approximations that we have intro- 
duced in the above algorithm, and it is of interest to ask 
how well the method works on simulated data. This also 
provides us an opportunity to demonstrate the advantages 
of this method for realistic observational scenarios where 
irregularly shaped fields complicate the calculation of the 
large-angle correlation function or power spectrum. These 
might arise from excising contaminated regions of the map 
or from sparse sampling strategies. 



4.1. Gauss 



Shear 



In Fig. ||, we show a Gaussian realization of a shear 
power spectrum that corresponds to the AC DM cosmol- 
ogy with source galaxies at z = 1 employed in White & Hr 
(2000). Gaussian distributed noise has been added to the 
625 ~ 8' x 8' pixels corresponding to 7™* = 0.4 and n = 56 
gal/arcmin 2 . The theoretical signal and noise power spec- 
tra are shown as lines in the left panel. We choose 7 ee, 
6 (3(3 and 6 e(3 bands according to the rules set down in 
the previous section. The exact /-ranges for the ee bands 
are given in Table U the other sets are similar save for the 
absence of the lowest /-band. Note that 1& c \a = 106 and 
Zpix = 2550. The recovered ee power spectrum is shown in 
Fig. H (top left) with sampling error estimates per realiza- 
tion obtained from averaging the Fisher variance estimates 
and run-to-run scatter over 200 realizations. Not shown 
are the (3(3 and e(3 bands where no statistically significant 
power was recovered. Notice that the two methods of es- 
timating the errors are in excellent agreement. The bands 
were chosen to be wide enough that their correlation due 
to survey geometry is negligible. For these Gaussian sim- 
ulations, the covariance matrix of the bands reflects this 
fact with negligible off-diagonal entries except for the end 
ee bands which have little intrinsic signal. 

4.2. N-body Shear 

In reality the shear field will be non-Gaussian due to the 
non-linearity of the underlying density fields. Because of 
the averaging effect of projection, the non-Gaussianity of 
the shear is much milder than that of the density field. To 
test its effects on the likelihood method, we use the simu- 



lations described in White & Hu (2000). We pixelize these 
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X (arcmin) 

Fig. 3. — Gaussian (top) and N-body (bottom) simulated shear (right) and recovered ee band powers (left) from 200 realizations. Solid lines 
represent the predicted power spectrum; dashed lines represent the assumed shot noise contribution. Errors are estimated from the Fisher 
matrix (shaded boxes) and the run-to-run scatter (open boxes). The simulated fields contain 625 pixels and are analyzed with 7 ee bands, 6 
j3/3 bands, and 6 e/3 bands. The latter two sets show no significant recovered power and have been omitted here for clarity. 
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Fig. 4. — Checkerboard sampled Gaussian shear and recovered power ee band powers from 100 realizations. The number and size of the 
pixels as well as the number and type of bands are the same as in Fig. § The mean power and errors are well recovered in spite of the more 
complicated sampling. 



simulations to the same level as used for the Gaussian runs 
and add the same amount of shot noise. A sample shear 
field is shown in Fig. || (bottom right). In the presence 
of this level of pixelization and noise, the non-Gaussianity 
of the N-body shear is in good agreement with th at esti- 
m ated fr om higher reso lution ./V-body simulations ( Jain et 



al. 20 DO; White fc Hu 2000) 



In Fig. | (bottom left), we show the power spectrum 
recovered from the simulated skies with error estimates as 



obtained in the Gaussian simulations. Non-Gaussianity 
does not bias the power spectrum estimates. Deviations 
in the lowest I bin may be attributed to finite box effects 
in the |Whitc fc Hr] (2000) simulations. 



Nevertheless the errors are slightly underestimated by 
the Fisher matrix in the intermediate regime where the 
intrinsic field is mildly non-Gaussian and the removal of 
shot noise does not dominate the errors. The increased 
variance mainly arises from the covariance of the Fourier 
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Table 1 



Band Covariance 



band 


37-90 


90-270 


270-450 


450-669 


669-1306 


1306-2550 


2550-5100 


37-90 


1.00 


-0.14 


0.23 


0.12 


0.13 


0.07 


-0.08 


90-270 


(-0.15) 


1.00 


0.08 


0.20 


0.17 


0.05 


0.02 


270-450 


(0.01) 


(-0.11) 


1.00 


0.37 


0.43 


0.25 


-0.15 


450-669 


(0.00) 


( 0.00) 


(-0.11) 


1.00 


0.38 


0.26 


-0.22 


669-1306 


(0.00) 


(-0.01) 


(0.00) 


(0.05) 


1.00 


0.50 


-0.38 


1306-2550 


(0.00) 


(-0.01) 


(0.00) 


(0.03) 


(-0.31) 


1.00 


-0.66 


2550-5100 


(0.00) 


( 0.01) 


(0.00) 


(-0.00) 


(-0.40) 


(-0.68) 


1.00 



NOTES. — Covariance matrix of the ee bands recovered from iV-body simulations. Upper: run-to-run covariance. Lower (parenthetical 
numbers): Fisher matrix estimates. The Fisher matrix underestimates the covariance in the intermediate regime where the signal is mildly 
non-Gaussian and dominates the shot noise. 
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Fig. 5. — Monitoring and removing excess noise. Top: Shot noise 
power underestimated by a factor of 100. The recovered 9 ee and 8 
j3P bandpowers (end 2 bands of each off scale) show an excess that is 
equal in the two bands and rises as I 2 to the pixel scale £ p j x f» 7500 
(see text). Bottom: Shot noise underestimated but removed with an 
added white noise parameter p no i se - Both the 7 ee bands (end band 
off scale) and the excess white noise, are well recovered. Boxes and 
shading represent errors as in Fig. H and each Gaussian realization 
of 60 utilizes 625 pixels. 



modes within the bands induced by non-linear mode cou- 
pling effects in the underlying density field. Mode cou- 
pling also correlates the bands themselves. We show the 
covariance matrix Gov a p/ (CovaaCoyg^) 1 / 2 of the ee band 
powers in Table [l] from the Fisher matrix and from the run- 
to-run covariance. Again in the intermediate regime, the 
covariance of the bands is underestimated by the Fisher 




matrix. 

Although not a severe effect, this underestimation sug- 
gests that in the absence of a sufficient number of fields 
from which the covariance matrix may be extracted di- 
rectly from the data, simulatio ns an d semi-analytic tech- 
niques (e.g. Scoccimarro et aL 1999; Cooray fc Hu 2000) 
should be used to calibrate the covariance matrix of band 
powers extracted from the data before using the measure- 
ments to constrain cosmological parameters. 

4.3. Irregular sampling 

One of the advantages of likelihood based methods is 
that they automatically account for any irregularity in the 
sampling or survey geometry, while maintaining an opti- 
mal weighting of the data on each angular scale. Sparse 
sampling techniques can be used to extend the dynamic 
range of power spe ctrum estimates for the same amount 
of observing time (Kaiser] 1998). Irregular sampling can 
also be used to test the effect suspect regions of the field 
might have on the results. 

For illustration we test the method on a large field sam- 
pled in a checker-board fashion with the same number and 
size of the pixels as above (see Fig. || right panel). Fig. || 
(left panel) shows the power spectrum recovered from 100 
realizations along with error estimates as above. Both the 
recovered power and the errors agree well. 

4.4. Excess Noise 

Including bands for the (3f3 and e/3 power spectra is 
not strictly necessary if all sources of noise have been ac- 
counted for in the noise covariance matrix. However they 
do provide a means of checking for any excess shot noise or 
systematic effects in the data. To illustrate this use, sup- 
pose that the initial estimate of the shot noise contribution 
were low by a factor of 10 (i.e. by 100 in power). As shown 
in Fig. [5] (top), the resulting power spectra in ee and (3 ft 
show excess contributions which scale as I 2 in bandpower 
and are equal in ee and f3(3 above the pixel scale. Since 
a white noise spectrum is not well approximated by a flat 
bandpower, errors in the pixel window function are exac- 
erbated leading to deviations near the pixel scale. 

In the event that excess noise is detected and that its 
correlation function or power spectrum may be parame- 
terized, the likelihood technique can be easily modified to 
include and effectively marginalize these noise parameters 
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from the data itself. The white noise case above provides 
the simplest example where there is only one extra param- 
eter pnoisc, providing an addition to the signal covariance 
matrix of the form 



Pn 



(y)(ob) 



N; 



-OijOab 



(24) 



In Fig. U (bottom), we show the result of dropping the 
(3(3 bands in favor of a white noise parameter. Both the 
true signal and excess white noise are well recovered by the 
technique. The parameter p no i se is plotted as a power spec- 
trum Cf olsc = p lw ise/n with error estimates from both the 
run-to-run scatter and the Fisher matrix. The covariance 
or degeneracy between the two is negligible since white 
noise has equal power in the two modes whereas the true 
signal does not. Excess noise from systematic errors will 
be more challenging to model, but the same principles and 
methods apply. 

5. CONCLUSIONS 

On large angular scales the shear field induced by weak 
gravitational lensing of background galaxies by large-scale 
structure is close to Gaussian. In this regime the relevant 
information is encoded in the angular power spectrum, 
l(l + l)Ci/(2ir). We have suggested a technique, commonly 
used in CMB analysis, for determining the angular power 
spectrum in this regime: the use of "band-powers" ex- 
tracted from the data by an iterated quadratic estimator 
( Bond, Jaffc fe Knox 1998) of the maximum likelihood so- 
lution. This technique has the advantage of automatically 
taking into account irregular survey geometries and vary- 
ing sampling densities. It provides an optimal estimate 
of the power spectrum which makes efficient use of all of 
the data on the relevant angular scales. We have tested 
the technique against simulated Gaussian and realistically 
non-Gaussian data, regular and irregularly sampled data, 
and with known and unknown amplitudes of shot noise 
from the intrinsic ellipticities of galaxies. In all cases, the 
mean band powers are recovered correctly. 

The technique introduced here is a result of one of many 
possible cross-fertilizations of CMB and weak lensing re- 
search. Indeed the experience gained in measuring the 
shear power spectra from noisy windowed data may feed 
back into the analogous problem for future CMB polar- 
ization studies. Techniques for handling large CMB data 
sets where the likelihood algorithm used here become s pro- 
hibitively expe nsive (e.g. Wandelt, Hi von fc Gorski 2000; 



Szapudi et al. 2000) will also be useful to lensing studies 



as the lensing fields become ever larger. 
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